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ABSTRACT 


This thesis presents a comparison of the measured and calculated 
flow field properties of a nozzle of unusual design which was used to 
produce an incompressible, low Reynolds number jet. The "nozzle" is 
essentially a porous metal plate which covers the end of a pipe. 
Results are presented for nozzle Reynolds numbers from 50 to 1000 
with velocities of 100 or 200 ft/sec. The nozzle produces a uniform 
velocity profile at nozzle Reynolds numbers well below those at which 
conventional contoured nozzles are completely filled with the boundary 
layer. A jet mixing analysis based on the boundary layer equations 
accurately predicted the flow field over the entire range of Reynolds 


numbers tested. 
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INTRODUCTION 


In the past decade several methods of measuring the temperature 
of the earth's upper atmosphere have been developed. In one of these 
methods a small sounding rocket is used to carry an instrument package 
and parachute aloft. The instrument package descends through the 
atmosphere suspended "beneath the parachute and telemeters temperature 
data to the ground. Due to the lov atmospheric density, the 
temperature transducer has a lov degree of convective coupling to the 
atmosphere. Solar radiation represents a significant portion of the 
total heat transfer vhich results in degraded accuracy. In addition 
the recovery factor must he knovn in order to correct for aerodynamic 
heating. Therefore interpretation of the data requires calibration 
on the ground with flov conditions equivalent to those experienced 
in the upper atmosphere. 

In a typical calibration arrangement, the temperature transducer 
is placed in a vind tunnel or jet of air to simulate the fall through 
the atmosphere. The flov Mach number and Reynolds number are dupli- 
cated to maintain the proper heat transfer characteristics. 

Hovever, for very lov Reynolds numbers, vhich are required to 
duplicate the conditions in the upper atmosphere, conventional vind 

tunnel nozzles develop very thick boundary layers. This problem 

/ 

becomes serious at nozzle Reynolds numbers of about 1000 and gets 
progressively vorse until the entire nozzle is filled with the 
boundary layer at Reynolds numbers of about 200 [l]. 
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Due to this Reynolds number limitation, a small research program 
v&s initiated to develop a low Reynolds number nozzle for a calibration 
facility at the Langley Research Center. A literature search revealed 
•that a multi-jet nozzle was a potential solution to the large boundary 
layer problem. In a multi- jet nozzle the flow is broken down into a 
number of sral.l jets which combine through viscous interaction to form 
a uniform flow downstream of the nozzle. A supersonic multi-jet nozzle 
(consisting of 37 holes in a 2 inch diameter steel plate) was tested at 
Ames Research Center with only partial success [l]. Although the nozzle 
did’ produce a uniform velocity profile, the flow was not established until 
it had gone many diameters downstream. 

The nozzle used in the present study is an extension of this 
concept to a large number of very smal 1 jets. With very small jets, 
local viscous interaction is much greater and a uniform flow is 
established almost immediately downstream of the nozzle. Since it 
did not seem feasible to accurately drill a large number of small 
holes, a number of naturally porous materials were tested on a trial 
and error basis. The first materials tested were paper filters with 
various porosities. All of these exhibited large deflections at 
fr-i gh flow rates which tended to distort the velocity profile. The 
addition of copper reinforcing screens produced adequate nozzle perform- 
ance except at very high flow rates. The nozzle material which finally 
evolved consists entirely of stainless steel screens which were 
sintered to form a 1/8 inch thick plate. This was very rigid 
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and produced a reasonably uniform velocity profile over a wide range 
of flowrates. 

The purpose of this thesis is to describe the performance 

i 

characteristics of this nozzle for nozzle Reynolds numbers between 
50 and 1000 and to compare the nozzle flowfield data with the calcu- 
lated results of an incompressible boundary layer type analysis. 

APPARATUS 

Figure 1 shows a schematic diagram of the test apparatus used 
in this investigation. The porous plate nozzle was bolted to the end 
of a 7-5/8 inch diameter nozzle pipe which extended through the wall 
of a large (55 feet diameter) vacuum chamber. The porous plate 
nozzle was made of 5 layers of stainless steel screen which were 
sintered to form a rigid plate 1/8 inch thick. The innermost screen 
had a mesh size of about 5 microinches. The outer screens were of a 
larger mesh size for strength. The sintered screens are sold commer- 
cially by the Eendix Corporation, Filter Division under the brand 
name Poroplate. 

Outside the vacuum chamber an air supply was passed through a 
dryer, pressure regulator, vertical tube flowmeter, and manual control 
valve into the nozzle pipe. The air stream stagnation temperature 
was measured inside the nozzle pipe using a shielded thermocouple. 
Downstream from the nozzle the flow Mach number was determined from 
measurements of the static pressure and the difference between the 
total and static pressures using the one-dimensional isentropic flow 
equations. 


— Vacuum cylinder vail 


Air 

supply 



Figure 1.- Schematic diagram of test apparatus. 
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The total pressure vas determined from impact tube measurements 
after applying a viscous correction. This correction is described in 

detail in the section entitled "Impact Tube Viscous Corrections." A 

j 

Baratron differential pressure transducer with a range zero to 10 
torr was used to measure the difference between the impact tube 
pressure and the static pressure. The static pressure throughout the 
jet was assumed to be the same as the background gas pressure since 
the flow vas subsonic. The background pressure vas measured at a 
point near the nozzle using a Baratron absolute pressure transducer 
with a range zero to 10 torr. Both pressure transducers were kept in 
a controlled temperature environment to minimize temperature effects. 

The flovfield surveys were made by moving the impact tube with 
a two-dimensional survey device. The survey device vas constructed 
and aligned so that surveys could be made either axially along the 
nozzle centerline or vertically along a nozzle radius. Calibrated 
potentiometers were used to indicate the distance in each, survey 
direction. 

All position, pressure, and temperature data were recorded on a 
Vidar digital data acquisition system. This system recorded at a 
rate of ^0 data channels' per second and produced either magnetic tape, 
printed output or both. 

/ IMPACT TUBE VISCOUS CORRECTIONS 

In this study the flow Mach number vas determined from the total 
end static pressures using the isentropic flow equations. The total 
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pressure vas obtained from impact tube measurements after applying 
viscous corrections. The impact tube used in this study is shown 
schematically in Figure 2. It vas constructed to duplicate (as nearly 
as possible) a probe described by Sherman [2]. In [2], several probes 
were calibrated to determine the magnitude of the errors resulting from 
viscous flow about the probes at low F.eynolds numbers. 

A discussion of impact tube errors at low Reynolds numbers and 
an analytical solution for certain types of probes is given by Schaaf 
[3]» For a probe pointing into the gas stream, a boundary layer forms 
at the stagnation point on the tip of the probe. Using the coordinate 
system shown in Figure 3, the y-conrponent of the Navier-Stokes 
equations for incompressible flow along the stagnation streamline can 
be written in the following form 





H 




(i) 


This equation can be integrated through the boundary layer to give 


P i = 


+ 


p u f/ 2 - 


(~) 

V 


= 6 


vl p 


y = o. 


( 2 ) 


since the pressure and velocity at the outer edge of the boundary layer 

/ 

are the same as in the free stream. From the continuity equation and 
flow symmetry about the stagnation streamline, it follows that 

3v _ 3u 3v 3u / v 

3y 3x - 3z ~ ” 3x 
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As a result of the boundary condition at the wall (^) 


equation (3) yields 


/dv» _ . 

( 3y y = 0 0 


Therefore equation (2) becomes 


* P H * 2 “ y=6 


From a~p-ote£tinl ^lovnenalysis 'the velocity gradient can 

be determined sjialytically for simple shapes , such as a spherical 
tipped probe [i;J. For a more complicated geometry such as the open- 
ended prcbe used in this study, ("|“~) g can approximated by 
k/4 (u^/D) where k is a constant which depends only on the probe 
geometry. The equation for impact pressure then has the form 


p p uj J 2 -f ^ y/2D 


fj_ — / v k U — I 

<°u l/z ~ 1 


Thus for large Reynolds numbers the correction term — is small and 

l\6 

p 

p. - p = pu /2. Vhen the Reynolds number becomes very small, 

2 

p^ - p^ becomes greater than pu ro /2. 
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As previously mentioned the constant, k, depends on the prohe 
geometry, and must he determined experimentally. Since the prohe used 
in this study was essentially a duplicate of one of the prohes tested 
by Sherman [2], a value of k equal to 6, determined hy fitting a 
curve through the data in Figure 5 of [2], was used. 

Once the value of k is known, the flow velocity may he 
calculated as follows for incompressible flov: 

R -R. -(/ +%)(PU*/Sj4+ eai ^y 

By rearranging and dividing hy p/2, the above equation becomes 

u a = Q 

from which one can solve for the velocity as 



The values of p and p in equation (9) are determined from the 
static temperature and pressure. Equation (9) therefore provides 
the corrected velocity in terms of the measured parameters. 

TEST CONDITIONS AND PROCEDURE 

Flovfield surveys were made at nozzle Reynolds numbers of 50, 
100, 200, 600, and 1000. Nominal flow velocities were either 100 


(7) 


( 8 ) 
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or 200 ft/sec. The resulting Mach numbers vere approximately O.O 89 
and 0.178, respectively, vbich vere sufficiently lov to insure essen- 
tially incompressible flov conditions. The Reynolds number was varied 

i 

by changing the density of the air. The flov conditions for each 
Reynolds number are listed in Table 1. 

t 

For testing at lov pressures a continuous flov capability is 
highly desirable in order to allov sufficient run time to establish 
stable flov conditions and make the required measurements. In the 
present study continuous flov conditions vere maintained for all 
nozzle flowrates . 

Before each test the static pressure transducer vas checked 
against a reference transducer to check for zero drift. The 
differential pressure transducer which vas used to measure the impact 
tube pressure vas checked for zero drift vhen the chamber vas at the 
desired static pressure just prior to establishing the flov. The 
desired flov conditions vere established based on the indicated dynamic 
pressure from the impact tube measurement vith precomputed viscous 
corrections. Preliminary surveys vere made to establish the fact that 
the jet vas axisymmetric. Surveys vere then made from the jet center- 
line outvard # along a radius at various positions along the jet axis. 
Survey data vere recorded cn the data acquisition system. After a 
complete flovfield survey, the process vas repeated for the next 
Reynolds number. 
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TABLE 1. TEST CONDITIONS 


ITozzle 

Reynolds 

Number 

Static 
Pressure 
x 103 cmHg 

Dynamic 
Pressure 
x lO^ mnHg 

Pi - P» 

3 

x 10 innHg 

Nominal 
Velocity 
ft /sec 

50 

100 

0.55 

1.86 

100 

100 

200 

1.09 

2.kl 

100 

200 

200 

k.k2 

7.06 

200 

6oo 

600 

13.26 

15-90 

200 

1000 

1000 

22.09 

2k. 13 

200 









JET FLOWFIELD ANALYSIS 


The Jet flovfield vas analyzed using the computer program 
described by Fox, Sinha and Weinberger [5]* The program has been 
adapted to the CSC 6000 series computer system and vas used vith only 
superficial modifications. 

Basically the program solved the compressible, axisymmetric 
boundary layer equations in the Von Kises plane using an implicit 
finite difference numerical technique. Thus the applicable equations 
for the f lev field shetehed in Figure 1 -are, as follows: 


Conservation of mass 


a C PU y) 

5 X 


*y 


- o 


Conservation of momentum 


( 10 ) 


_ cj M 3 LA 


-1l -t±i 

ax y *y 
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(n) 


Conservation of energy 


. j i 
c r d x 


fVc, 




il 


- ns? 

U 3X + 


« w 

/ *y 


y) + Cc02. 9>. ^ y [ 3 / J 


( 12 ) 
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Conservation of species 


0 _ 2 S; 

r t / J7 - r““ v 


/c yj± f<£ js, 

- C * 0 ^ ^ y <iy ( 13 ) 

c ~ 1,2 

Initially, the velocity, temperature, and gas concentrations are 
specified as: 


U [*t> j y) 
d"/ (X&jf) 


- U a (?] 

~ T* (y) 


( 14 ) 


The boundary conditions at the center of the jet due to flow symmetry 


are 


^ u (x , c) - 6^£> - ~ Q (15) 

*y <) / 


the cut 

er edge cf the jet. 

J im 

U (*> y) - 'Je. 

y- 9- 00 

J * Hi 

T (*■,/) * T*. 

y-*? «■=' 


^ ^ 



( 16 ) 


By defining dinens ionless quantities as 

X - x/R 
y = y/R 
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a = u /tv t 
v = v/ u t 

e - f / q 
f - r/r c 

9 * 


(17) 


t&e conservation equations can be solved in the modified Von 2dises 
plane vnich is defined for axi symmetric flov by: 




r e/X 

j ^ 


( 18 ) 


and 



(19) 


It should be noted that since the edge conditions were used to 

normalize the velocities, the case of u = 0 cannot be solved. 

e 

However, this is usually not a problem since u may be made very 

e 

small, and was not zero in the actual experiment. 


DATA ACCURACY 


An error estimate vas cade for each of the variables of interest 
in this study. The errors in the measured quantities vere determined 
directly free instrument characteristics. The errors in calculated 
quantities vere based on the errors of each measured quantity which 
vas used in the calculation. An error summary is presented in Table 2. 

The Eass flcv and total teeperature accuracies are based on the 
manufacturers' specifications for each instrument. Since the nozzle 
Eass flcv rate a?as uised.^Rly for a .qualitative check of the average 
nozzle velocity, an additional mass flow calibration did not appear to 
be justified. The +2°? total temperature error represents only 0.^ 
percent of the absolute temperature and is a relatively small error 
source. 

The error in hoth the static pressure and the difference between 
the impact tube and static pressure is primarily due to instrument 
zero drift. The static pressure transducer was located in a small 
controlled environment chamber inside the large vacuum chamber. The 
transducer vas not readily accessible for zero calibration since this 
required physically pumping the transducer to zero pressure. However, 
a reference transducer which vas calibrated against a secondary 

standard vas kept pumped down to essentially zero pressure. Just 

/ 

prior to each test, the static pressure was measured using hoth 
transducers end the reading of the static pressure transducer corrected 
to the value indicated by the reference transducer. Zero check for 


IT 
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TABLE 2. DATA ACCURACY 



Re 


+2? to +11? 


Based on errors in p, v, and y 
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the differential pressure transducer vas relatively easy. This was 
-done just before each test before flow conditions were established. 

The error in measuring the survey probe position was estimated 
to be about +. 03 inch primarily due to gear backlash on the survey 
device. The error was estimated by clamping the probe in a fixed 
position and monitoring the position readout while turning the motor 
in either direction. 

Table 2 also gives an error estimate for variables calculated 
frcm the primary measurements. The density error is based on a 
combination of the temperature and pressure errors through the 
perfect gas equation assuming that the gas constant is known. The 
viscosity error is based on the temperature error and the Sutherland 
viscosity equation given in [6]. 

The error in determining the flow velocity is based on a 10 
percent uncertainty in the impact tube viscous correction constant, 
the temperature error for determining the speed of sound, and the 
pressure errors for determining the Mach number. The error in the 
speed of sound is only +0. 2 percent assuming that the ratio of 
specific heats and gas constant are known. The Mach number error in 
terns of the pressure ratio error is given by 

hjfjrjn rf 

J (£) ( 2o) 



or in terms of the measured quantities 
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The velocity error is then the sun cf the Mach number error and speed 
of sound error. 

The error in Reynolds number is sinply the sun of the errors in 
density, viscosity, and velocity. The error in measuring probe 
diameter is included in the uncertainty of the impact tube viscous 
correction constant. 

RESULT'S PSD DISCUSSION 

Figure 5 Tdraws the ’davsicpsaent of the Jet velocity profiles in 
the dovn stream direction for a nozzle Reynolds number of 1000. 
Velocity profiles are shown for three axial stations: x/R = 1, 2 . 5 > 

and 7- For each velocity profile, the ratio cf local velocity to 
the maximum jet velocity, v/ v^^ , vas plotted against the non- 
dimensional radial coordinate, y/R. The symbols represent the 
measured data and the solid line represents the computed results . 

The measured and computed results are in good agreement, especially 
in the mixing region except at x/R = 7- As can be seen a very 
large potential core exists as far downstream as x/R = 7 which 
vas the maximum distance for which measurements were taken. The cal- 
culated results indicate that the potential cere would extend to about 
x/R = 30. 

Figure 6 shows a set of velocity profiles for a nozzle Reynolds 
number of 6C0. Again, v/v_ is plotted against y/R for x/R = 


1, 2.5, and 7.79- The agreement between calculated and measured 



Figure 5.- Nozzle velocity profiles for Re = 1000 
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results is good except at 7.79. The potential core is still reasonably 
large at x/R = 7.7 9 and according to the analysis, it extended down- 
stream to about x/R = 20. 

Figure 7 shows typical velocity profiles for a nozzle Reynolds 
number of 200. v/v is plotted against y/R for x/R = 0.5, 2.5, 
and 6.0. The agreement between measured and calculated results is 
again very good. The potential core extends downstream to about 
x/R =6. At x/R = 0.5 and 2.5, there is a reasonably large region of 
uniform velocity. It should be noted that at Reynolds numbers of about 
200, a conventional 'contoured nozzle would be nearly filled with 
boundary layer at the nozzle exit. 

Figure 8 shews the calculated and measured flew field data for a 
nozzle Reynolds number of ICO. v/v ^ is plotted against y/R for 
x/P. = 0.5, 1.0, and 2.5. The potential core extends downstream just 
slightly beyond x/R = 2.5. At x/R = 0.5 and 1.0 there still remains 
a relatively large area of uniform flow. 

Figure 9 shows the velocity profiles for a nozzle Reynolds number 

of 50, the lowest Reynolds number at which flow surveys were made. 

v/v is again plotted versus y/P. for x/R = 1.0, 2.5, and 5«0. 
nax 

The potential cere extends downstream to only about x/R = 1.0. The 
calcinated results, which agree with the data very well further .♦ 

downstream, indicate that there is a large region of uniform flow 
Just downstream of the nozzle. For example, at x/R = 0.5, the 
calculated results indicate that the potential core extends out 
radially to y/R = 0.5. 
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As mentioned, earlier, the agreement between the measured and 
calculated results is very good with only a few exceptions. The 
agreement would probably be even better if it were not for two prob- 
lems in calculating the flow field. The first problem is with the 
initial velocity profile which is used to start the flow field 
calculations. It was assumed that at the nozzle exit the velocity 
was uniform across the entire nozzle. At the edge of the nozzle, 
it was assumed that the velocity made a step change to the value 
associated with the background gas. However, the data indicate that 
the velocity was not quite uniform across the nozzle, especially at 
the higher Reynolds numbers. These errors enter the calculations 
and produce a large part of the discrepancies. 

The other problem is in the value of velocity used for the 
background gas. The computer program which was used for these 
• calculations is designed to handle either the flow field of a Jet 
or the wake flow behind a body. However, wake flow calculations are 
its primary function. In wake calculations, it is very logical to 
nondimens ionalize the velocity in the wake to the external flow 
velocity, u . However, this precludes calculating the flow field 
of a jet flowing into a fluid at rest since u g would be zero and 
the nondimensional velocities would be infinite. This normally would 
not be a problem since u g may be made arbitrarily small, and, as Pai 
[7] has shown, the solution is not overly sensitive to changes in u . 

1 ® 

However, for very low nozzle Reynolds numbers there is an additional 
problem. The numerical calculation step size is based on the Reynolds 
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number of the external flow. At low nozzle Reynolds numbers, small 

: values of u produce prohibitively long run times. In this 
e 

study, u g was taken to be 5 percent of the jet velocity for all 
calculations . This produced reasonable run times without introducing 

■ excessively large errors in the calculated velocity profiles. 

j 

| The primary effect of u is to change the jet spreading rate. 

I 6 

I 

j This effect is shown in Figures 10 through l4. The 0.5 velocity 

! radius (value of y/R at which v/v = 0.5) is plotted versus x/R 
j max 

I 

' for Reynolds numbers of 1000, 600, 200, 100, and 50. The 0.5 velocity 
radius is a measure of the jet width and its change in the downstream 
direction is a measure of the jet spreading rate. ■ These figures show 
that the experimentally determined spreading rate is greater than that 
| -calculated using a value of u g which is 5 percent of the jet velocity. 

At nozzle Reynolds numbers of 50 and 100, the nozzle velocity was 
100 ft/sec rather than 200 ft/sec which was used for all the other 
; tests conditions. In order to determine if this change had any 

i affect on the results, flow field calculations were made for both 

i 

! 100 ft/sec and 200 ft/sec velocities at the same nozzle Reynolds 
j number. The nondimens ional velocity profiles were identical which 
! indicates that the mixing and spreading of an incompressible jet is 

j 

a function of only the Reynolds number . 

j * 

i 

» 


j This is also what one would expect intuitively. 



velocity radius for Re * 1000 
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Re = 100 
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CONCLUSIONS 


1. The porous plate nozzle produced a reasonably uniform 

velocity profile over a range of Reynolds numbers from 1000 down to 
50. This contrasts with a conventional contoured nozzle which would 
have a large boundary layer at a Reynolds number of 1000 and would be 
completely filled with boundary layer at Reynolds numbers on the 
order of 200. . 

2. A conventional, boundary layer type analysis was sufficient 
to accurately calculate the Jet flow field for nozzle Reynolds 
numbers as low as 50. 

3. The calculated mixing and spreading of an incompressible Jet 
issuing into a medium at rest was a function of the Reynolds number 
only. Calculations at the same nozzle Reynolds number but different 
velocities produced essentially identical results. 
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